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We study double-exchange models with itinerant tig electrons in spinel and pyrochlore crystals. 
In both cases the localized spins form a network of corner-sharing tetrahedra. We show that the 
strong directional dependence of t^ g orbitals leads to unusual Fermi surfaces that induce spin super- 
structures and noncoplanar orderings for a weak coupling between itinerant electrons and localized 
spins. Implications of our results to ZnY^CU and Cd20s207 are also discussed. 



Orbital degrees of freedom have attracted much atten- 
tion due to their crucial role in the stability of many 
unusual phases of correlated materials In particular, 
the presence of degenerate orbitals in frustrated magnets 
can lift the spin degeneracy through various spin-orbital 
interactions. For Mott insulators with spins residing on 
a frustrated lattice, such as kagome or pyrochlore, geo- 
metrical constraints prevent spins from reaching a simple 
Neel order. The occurrence of long-range orbital order 
due to either Jahn- Teller distortion or orbital exchange 
reduces the magnetic frustration by creating disparities 
between nearest-neighbor (NN) exchange constants and 
paves the way for magnetic ordering. 

However, some of the magnetic orders observed in ge- 
ometrically frustrated compounds are difficult to explain 
starting from the strongly coupled Mott- insulator regime. 
For example, several vanadium spinels exhibit a 

complicated magnetic structure with ttii * ' * collinear 
ordering along certain chains that is very puzzling from 
the viewpoint of localized spin models. Below we shall 
provide a simple explanation for the observed spin super- 
structures based on a double-exchange (DE) model which 
takes into account orbital degeneracy. The DE model 
arises naturally for multiband compounds in which a nar- 
row band of localized electrons coexists with a wider band 
of itinerant electrons. It can also be viewed as a mean- 
field approximation to the Hubbard Hamiltonian. A well 
studied case is the DE model with itinerant e g electrons 
on the cubic lattice [H|. This model has been shown to 
describe the rich physics of colossal magnetoresistance in 
perovskite manganites. 

Recently there has been tremendous interest in DE 
models on frustrated lattices [6-131. The Fermi surface 
geometry plays a crucial role [8|, LL3| in the nonlocal ef- 
fective spin-spin interaction that results from integrating 
out the itinerant electrons in the weak-coupling regime. 
The magnetic structures stabilized by itinerant electrons 
are thus often difficult to understand using short-range 
spin models. For example, an unusual noncoplanar mag- 
netic order, in which spins on different sublattices point 
toward the corners of a tetrahedron, is shown to appear 
in different coupling regimes and various commensurate 
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FIG. 1: A unit cell of the pyrochlore lattice and the config- 
uration of local oxygen octahedra in (a) spinel and (b) py- 
rochlore crystals. The blue and red balls denote the (B-site) 
transition-metal and oxygen ions, respectively. 



filling fractions on the triangular lattice (j- 11|. Recent 
investigations of DE models on pyrochlore lattice also re- 
veal interesting behaviors such as electronic phase separa- 
tion 12[ and a complex noncoplanar order [13| at quarter- 



filling. However, most of these studies ignore the orbital 
dependence and consider only isotropic electron hopping. 

In this Letter we examine DE models with itinerant 
t2 g electrons in both spinel and pyrochlore structures 
with general formulas AB2O4 and A2B2O7, respectively. 
We consider localized classical spins residing on the B- 
sites of the crystal which form a 3D network of corner- 
sharing tetrahedra. The 06 octahedron surrounding the 
B-sites creates a cubic crystal field which splits the d- 
levels into the e g doublet and the lower-energy t2 g triplet. 
The strong dependence of electron hopping on orbital ori- 
entation leads to peculiar Fermi surfaces in both cases. 
In particular, the electron subsystem reduces to a set of 
cross- linking Kondo chains in spinels. We show that this 
feature leads to a weak-coupling instability towards the 
previously mentioned ttii superstructure in vanadium 
spinels. Fermi surface nesting of different origin leads to 
noncoplanar all-in-all-out magnetic order in pyrochlores, 
which is a candidate state for the intermediate insulating 
phase of Cd20s207. 

Frustrated Kondo- chains in spinels. In the spinel 
structure, a common quantization axis can be defined for 
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t2 g electrons at all crystal B-sites [Fig. 0Ja)]. The shape 
of the t2 g orbitals is such that the strongest overlap is 
between the same orbitals along a particular NN direc- 
tion, e.g., between two d xy orbitals along either a [110] 
or [110] bonds in the xy plane. Keeping only this domi- 
nant term, electrons in a given orbital state can only hop 
along the corresponding (110) chain on the pyrochlore 
lattice. We thus divide these chains into three types, 
yz, zx, and xy, depending on the active orbitals along 
the chain. Since the kinetic energy preserves the orbital 
flavor, the Hamiltonian is a sum of contributions from 
different orbital sectors: H = H m , where 

H m = ~t (4ma C jma + h ' C ') " Jh ^ S * ' S *™> W 

(ij)\\m i 

m = xy, yz, zx and s im = \ £ a/3 4 ma (T a ^c im ^ Here the 
first term describes NN hopping of t<i g electrons along a 
(110) chain of type m; t is the dominating ddcr transfer 
integral. cj~ mQ , is the creation operator for d-electrons 
at site i with orbital flavor m = xy,yz,zx and spin 
a =t, I- The second term in Eq. (pQ) describes an effec- 
tive on-site Hund's coupling between t 2g electrons and 
localized classical spins (with normalization |S^| = 1). 
Regarding model (pQ) as a mean-field approximation for a 
three-band Hubbard Hamiltonian that has the same ki- 
netic energy term as H, the effective coupling constant is 
J H = 4([//9 + 4J#/9)|(si)| where U + J H is the Coulomb 
repulsion between two electrons in the same orbital and 
Jh is the bare Hund's coupling [lij . 

H models a collection of ferromagnetic (FM) Kondo 
chains coupled together by the local moments. While a 
classical Kondo chain is a relatively simple system, the 
fact that each spin is shared by three chains with differ- 
ent orbitals introduces geometric frustration. Numerical 
methods such as Monte Carlo (MC) become very inef- 
ficient for conventional 3D DE models because the di- 
mension of the electron Hamiltonian to be diagonalized 
for each spin update scales as L 3 x L 3 for systems with 
linear size L. On the contrary, for H, one only needs 
to diagonalize matrices whose dimension scales as L x L 
for the three Kondo chains intersecting at the updated 
spin. The reduced dimensionality of the problem thus 
allows for studying the ground states of H with the aid 
of large-scale MC simulations. 

We first consider the case with three (i-electrons per 
site. The electron energy is minimized by placing one 
electron at each of the three different ID bands, giving 
rise to half-filled Kondo chains with a Fermi wavevec- 
tor kp = tt/21 [Fig. [2(c)], where / is the NN distance. 
The two Fermi points are nested by a commensurate 
wavevector qi/ 2 = 2kp = tt/1, leading to magnetic Neel 
order in the presence of Hund's coupling. However, di- 
rect inspection shows that such a collinear Neel order 
cannot be simultaneously attained on all chains of the 
pyrochlore lattice. Instead, MC simulations on L = 8 lat- 




FIG. 2: (a) The noncoplanar all-in all-out magnetic order 
and (b) collinear spin superstructure with wavevector Q = 
(0, 0, 2tv/cl) in spinels. Note the ttii ■ ■ ■ order along yz and 
xz chains and a titi ' • ' Neel order on xy chains in (b). Panels 
(c) and (d) show the corresponding electron band structures 
of the three ID chains with different orbitals. 

tices (with 16L 3 spins) show that the total energy is min- 
imized by the noncoplanar all-in-all-out spin order shown 
in Fig.(2fa). The magnetic order of each chain consists of 
FM and staggered components, which are perpendicular 
to each other. It is worth noting that any global rotation 
of this noncoplanar spin ordering leads to another ground 
state due to the SU(2) invariance of H. 

The situation is more complicated for transition met- 
als with two d-electrons per site like the vanadium spinels 
AV2O4 where A = Zn, Cd, or Mg. In ideal cubic spinel, 
equal distribution of electrons among the Kondo chains 
corresponds to 1/3 filling fraction. The classical ground 
state of a single Kondo chain at 1 / 3-filling has a tti 
magnetic order with a period of 3 1. Again such a simple 
arrangement of spins is precluded by geometric frustra- 
tion. Our numerical minimization on large finite systems 
yields a 3D noncoplanar magnetic order with wavevector 
Q = ^ Q, |, l), where a = 2^/21 is the length of a con- 
ventional cubic unit cell; the extended magnetic unit cell 
contains 108 spins. 

Collinear superstructure in vanadium spinels. Instead 
of the above complex order which preserves orbital de- 
generacy, experiments showed that vanadium spinels un- 
dergo a cubic-to-tetragonal structural transition with lat- 
tice constants c < a = b [l-Q- Contrary to the elon- 
gation, which is favored by a Jahn- Teller ion with two 
tig electrons, the observed tetragonal compression can 
be understood as originating from the band Jahn- Teller 
instability. The lattice distortion results in a crystal- 
field splitting of the t<i g levels as shown schematically in 
Fig.Efd). With two <i-electrons per site, the lower energy 
xy orbital is always occupied by one electron. Since the 
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lattice distortion also reduces the hopping integral of xy 
electrons, we assume that the d-electrons in the xy or- 
bital are localized and constitute the local moments {S^}. 
Therefore, we consider the DE Hamiltonian: 



DE 



^ Hm + JAF ^ Si-Sj. 



(2) 



m=yz,zx 



(ij)\\xy 



Here J af — &t y /U is the exchange constant between 
localized spins along the [110] and [110] chains. As the 
magnet is cooled, antiferromagnetic (AFM) spin correla- 
tions develop first along these chains as indeed observed 
[H. A long-range 3D magnetic order resulting from in- 
teractions between different spin chains sets in at a lower 
temperature [2-4]. However, the crossing-chain coupling 
is geometrically frustrated if only NN spin interactions 
are taken into account [TBI, [l6[ . 

Here we provide a simple picture of the unusual mag- 
netic order of these vanadates based on the DE model (|2]). 
Because the other d electron can occupy either yz or zx 
orbitals, the corresponding bands are both 1/4 filled. In 
presence of Hund's coupling, J#, the usual Fermi-point 
nesting thus leads to the formation of ttll ' ' * super- 
structure with qiu = 2k f = tt/21 on both yz and zx 
chains. In the weak-coupling regime, this collinear or- 
dering (same amplitude for ±(71/4) is always more sta- 
ble than the single-g spiral order because both wavevec- 
tors, qi/4 and —#1/4, are required to gap the two Fermi 
points of each chain. The corresponding 3D collinear 
magnetic order [Fig. Efb)] characterized by wavevector 
Q = (0, 0, 1) is consistent with the experiments (we shall 
from now on express the wavevectors in unit of 2ir / a for 
convenience). The mechanism for the formation of this 
magnetic order is similar to the orbitally induced Peierls 
instability in the spinel MgTi204 [171 ]. 

Although the above collinear spin order can also be 
explained within a local spin picture, an ad hoc third- 
neighbor AFM exchange has to be introduced in order 
to stabilize the t1M4 structure along yz and zx chains 
[l5| . On the other hand, our approach based on the itin- 
erant DE model provides a natural explanation for the 
formation of these superstructures in absence of orbital 
order. In addition, recent ah initio calculations and ex- 
perimental studies indicated that some vanadium com- 
pounds are indeed close to the metal-insulator transition 
jlil . [HI, giving further support to the itinerant picture 
adopted here. Although the above conclusion is valid 
only for weak J#/£, and is based on a mean- field treat- 
ment of the multi-band Hubbard model, a more exact 
calculation that takes into account the electron correla- 
tions gives a consistent result which will be presented 
elsewhere. 

Noncoplanar magnetic order in metallic pyrochlore. 
We now turn to DE model with degenerate orbitals on 
the pyrochlore structure [Fig. Gib)]. Our theory provides 
a plausible explanation for magnetic ordering and metal- 
insulator transition in the pyrochlore oxide Cd20s207. 



This compound undergoes a continuous metal-insulator 
transition at Tmi ~ 225 K [20|422| - The resistivity in- 
creases by 3 orders of magnitude upon cooling below 
Tmf The transition is accompanied by a sharp reduc- 
tion of magnetic susceptibility, indicating the occurrence 
of AFM order [2l|. The specific- heat anomaly at Tmi 



is found to be well described by a mean-field BCS-type 
phase transition. The electron activation energy obtained 
from resistivity measurements also exhibits a BCS-like 
behavior near Tmi [H|. 

These experimental observations justify a mean- field 
approach for the metal-insulator and magnetic transi- 
tion in Cd20s207. As discussed above, the mean- field 
approximation reduces the multiband Hubbard model to 
the following DE model: 



H 



MF 



ran t r 

ij ima jna 



JH ^ ^ Si • Si; 



.■(3) 



Here the orbital index m refers to the quantization axes 
of the local crystal fields which are different in the four 
nonequivalent crystal B-sites [Fig.QJb)]. Contrary to the 
case of spinels, the orbital flavor is not conserved by the 
kinetic term. To obtain the hopping matrix, we expand 
the t2 g orbital wavefunction at a given sublattice s in the 
basis of common coordinates for the cubic pyrochlore: 
\<t>m) — a m/el^)- The details of the transformation coef- 
ficients a s mk can be found in Ref. [23j . The resulting hop- 
ping matrix is t™? = Y.ki a ™k a n\ (M H Mi) • Here the 
transfer integral ((t)k\H t \(t)j) is expressed using the Slater- 
Koster (SK) parameters [24[. 

We again start by considering only the dominant dda 
hopping in the SK parameters; the calculated tight- 
binding spectrum is shown in Fig. [3^a). In the metallic 
pyrochlore Cd20s207, the Os 5+ ion has three d electrons 
corresponding to a half-filled band. The Fermi level lies 
at e = for this filling fraction and the resultant Fermi 
'surface' consists of three lines and four points at the 
boundary of the Brillouin zone. The three Fermi lines 
are diagonals of the square surface at the zone boundary, 
while the Fermi points are located at the high symmetry 
L-point k L = (i, \, |) [see Figs. [3^ a) and (b)]. 

Interestingly, this unusual Fermi surface can be nested 
by three wavevectors Qi = (1,0,0), Q2 = (0,1,0) and 
Q3 = (0,0,1). In particular, the Fermi lines are topo- 
logically equivalent to three 'circles' each of which can 
be completely nested by one of the Q vectors (l3| . To 
determine the optimal ground state we minimize the en- 
ergy among all the spin orderings for which the non- 
interacting system has a divergent susceptibility. In other 
words, we introduce a variational amplitude for the uni- 
form ordering with Qo = and each wavevector 
that leads to perfect nesting of the Fermi surface. Re- 
stricted to this particular set of magnetic structures, our 
simulated-annealing minimization yields a noncoplanar 
spin order characterized by a single wavevector Q = 0; 
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FIG. 3: (a) Band structure of the orbital-dependent tight- 
binding model in the pyrochlore crystal structure. At half- 
filling, the Fermi level is at e = 0. (b) First Brillouin zone of 
the fee lattice which is the underlying Bravais lattice of the 
pyrochlore crystal. The Fermi 'surface' for a half-filled band 
include three sets of Fermi lines (the XW segment in (a)) and 
Fermi points at the L points, (c) Band structure in the pres- 
ence of nonzero Hund's constant Jh — 0.05£. (d) Calculated 
density of states for zero and nonzero Hund's coupling. 

the magnetic unit cell is the same as the crystal one. 
Spins on the four inequivalent sites point toward the cor- 
ners of a tetrahedron. The so-called all-in-all-out struc- 
ture shown in Fig. [2fa) is an example of the noncoplanar 
'tetrahedral' order. The corresponding band structure is 
shown in Fig. [3^c) for Hund's coupling Jh = 0.05t. A 
charge gap opens at the original Fermi energy, as can also 
be seen in the calculated density of states [Fig. E(d)]. 

The noncoplanar spin structure obtained above can 
be a strong candidate for the magnetic order below 
Tmi in Cd20s207. This simple q = order also pre- 
serves the cubic symmetry. Experimentally, the metal- 
insulator transition was found to be accompanied by a 
slight change in unit-cell volume of less than 0.05%. More 
importantly, no change in crystal symmetry was observed 
below Tmi • Although the exact magnetic structure is yet 
unclear, the q = noncoplanar order is consistent with 
a recent /iSR measurement j22|. Interestingly, upon fur- 
ther cooling, an incommensurate spin density wave dis- 
continuously develops below T ^ 150 K j22j. This might 
indicate the breakdown of mean- field approximation deep 
in the insulating phase where strong electron correlations 
play a predominant role. 



In summary, we have studied the DE model with ti Q 
electrons on the pyrochlore lattice. By taking into ac- 
count the orbital-dependent hopping, we showed that 
magnetic properties of spinels close to the metal-insulator 
transition can be understood using the picture of cross- 
linking Kondo chains coupled by localized moments. Our 
theory provides simple and elegant explanations for the 
unusual spin superstructure observed in several vana- 
dium spinels. We also proposed a novel noncoplanar 
'tetrahedral' order for the magnetic insulating phase of 
the pyrochlore Cd20s207. 
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